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Abstract. Free quadratic harness is a Markov process from the class of qua- 
dratic harnesses, i.e. processes with linear regressions and quadratic condi- 
tional variances. The process has recently been constructed for a restricted 
range of parameters in [7] using Askey- Wilson polynomials. Here we provide 
a self-contained construction of the free quadratic harness for all values of 
parameters. 



1. Introduction 

Quadratic harnesses were introduced in [3] as the square-integrable stochastic 
processes on [0, oo) such that for all t, s > 

(1.1) E[X t ] = 0, E[X t X s ] = min(i, s), 

conditional expectations E [Xf \ J- S . u ] are linear functions of X s and X u , and second 
conditional moments E [Xf\ F StU \ are quadratic functions of X s and X u 

(1.2) E[Xf\T s , u ] =Qt, s ,u(X s ,X u ), 
where 

(1.3) Qt.s,u(x, y) = A t ^. u x 2 + B ti8jU xy + C t .s, u y 2 + D t<s<u x + E tjS<u y + F MiU , 

and At lS ,u, ■ ■ ■ j F t .s^u are some deterministic functions of < s < t < u. Here, J- S . u 
is the two-sided cr-algebra generated by {X r : r 6 (0, s] U [u, oo)}. We will also use 
the one sided er-algebras Tt generated by {X r : r < t}. 
It follows that for all s < t < u 

(1.4) E[X t \F s>u } = at, s , u X s + b t , s , u X u , 

with at jS ,u = (u — t)/(u — s) and bt, s ,u = (t — s)/(u — s), and, under certain technical 
assumptions, there exist five parameters 77, 9 £ K, a, t > and 7 G [—1, 1 + 
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such that 



(1.5) 



and 

(1.6) Var[X t |JF S)U ] 
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see [3 Theorem 2.2]. 

Quadratic harnesses may have orthogonal martingale polynomials (see [3] for the 
assumptions and |10| for some exceptions), some explicit examples of which have 
been worked out in Section 4 of [3] , for some of them the corresponding quadratic 
harnesses were constructed in a series of papers [1] , [S] , [5] . A recent development in 
proving the existence of quadratic harnesses is [7] , where the machinery of Askey- 
Wilson polynomials have been used to construct the processes for a wide range of 
parameters 77, 8, a, r and 7. 

However, in some cases, the theory developed in [7], brings some unnecessary 
limitations for the values of parameters assuring the existence of the given quadratic 
harness. One of them is the case of free quadratic harness defined in section 4.1 of 
[3]. Free harnesses have parameter 

(1.7) 7 = -ctt; 

their orthogonal martingale polynomials were identified in [3l Proposition 4.3], for 
a, t > 0, err < 1 and 1 + a(3 > 0, where 

T) + 9a j]t + 8 



(1.8) 



1 — ar ' 1 — err 

On the other hand, Corollary 5.3 in [7J, which discusses the range of parameters 
that guarantee the existence of the free harness, requires additional assumption 

(1.9) 2 + r/9 + 2<7T > 

to be compliant with the theory developed in [7J; nevertheless, the univariate Askey- 
Wilson distributions are still well defined when 2 + t]8 + 2ar < 0. 

The goal of this paper is to show that free quadratic harness exists without 
assumption ()1.9|) (this is stated in the main result of the paper - Theorem 1 1.1|) . The 
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technique we use is similar to the one used in previously mentioned work ([4], [5], 
[6]), although various details differ. We rely on explicit three step recurrences for the 
orthogonal martingale polynomials and on explicit connection coefficients between 
related families of orthogonal polynomials. We also use an operator representation 
to prove the quadratic harness property of the constructed process. The paper is 
self-contained and does not use any results from [7]. 
Our main result is 

Theorem 1.1. For a, r > 7 err < 1, 7 = —err, and 1 + a/3 > 0, there exists 
a Markov process (^t)te[o,oo) such that (iLl) . (fl~4"| and dL6l) hold. The process 
(X t )t£[o.oo) is unique among the processes with infinitely- supported one- dimensional 
distributions that have moments of all orders and satisfy (|1.4[) and (|1.6[) with 

the same parameters 77, 6, a, t and 7 = — ar. 

The proof of this theorem is given in Section |B] after all auxiliary technical results 
are established. 

It is worth mentioning that for some values of parameters, the univariate laws 
of free quadratic harnesses are the first component of a two-state free convolution 
semigroup (for the free bi-Poisson process case see [6] ; for an extension to the case 
of a = see [U Proposition 5]). 

2. Orthogonal martingale polynomials and (j-commutation equation 

This section presents a heuristic principle that can be used to find the recurrence 
for the polynomials orthogonal with respect to the conditional law £(X t \X s ). (For 
the proof of Theorem 11.11 such a derivation is not needed, as the actual recurrence 
used in the proof can be accepted as a guess.) 

In [3] we defined the orthogonal martingale polynomials associated with the 
process (X t ) t as martingale polynomials (i.e. such that 

(2.1) E[p n {X t ;t)\F s ] =p n (X 8 ;s), 

holds, whenever < s < t) that are orthogonal with respect to the one-dimensional 
distributions of the process. Theorem 2.3 from [5] states that their Jacobi matrix 
is linear in t, 



C t = 



( 7oi + <5 

Ott\t + Pi 
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oatf + p 2 








£2*4 
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and for n > 1 the coefficients a n , /3„, 7„, 5 n , 4> n , £n satisfy 
(2.2) a 2 ra n a n+ i + era nj 5„ + i7 + a(3 n /3 n+ i = aa n+ if> n 



(2.3) /3 



aa n+1 (i n +7„ + i)r+ (aa n+1 S n+ i + fi n +\'1n)l+ 



(2.4) /3„ + i£„ + i + 7„<5„ + <ra n ip n = (aa n+ i£ n+ i + 7^ + aa n e n )r+ 

(aa n+ iip n+ i + i n S n + /3 n £nh + Wn+lfn+1 + 8* + /3 n ip n )a + j n 6 + 5 n rj + 1, 
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(2.6) e n (p n+ i = e n e n+ ir + e n+ itp n ~/ + <p n tp n+ ia, 
with the initial values given by 

(2.7) ai =0, fc = 1, 70 = ^0 = 0, e x = 1, fa = 0. 
Now, for a fixed r > 0, conditionally on X r , the process (Y( )t>o 



- V 1+7^+^ {Xr+t ~ 

still is a quadratic harness. Therefore when one considers the conditional distri- 
bution C(X t \X r ) and the corresponding orthogonal polynomials Q n (y;x,t,r), [3j 
Theorem 2.3] implies that their Jacobi matrix is again linear in t and that its entries 
again satisfy relations (|2.2I) - (12.6[) . The only difference is that the initial values for 
the sequences should be modified as follows: 
(2.8) 

l + rjx + <7x 2 -T (l + i]x + ax 2 ) 

ax = 0, Pi = 1, 70 = 0, d = x, e\ = — , (pi = — — , 

1 + ar 1 + err 

as we choose the first polynomials Q n as 

Q-x(y;x,t, s) = 0, Qo(y;x,t,s) = 1, Qi{y, x, t, s) = y - x, 

and 

Q-z{y,x,t,s) = y 



1 + at 

f (l + as 

y< x 



(1 + as)(l + 7) (1 + as) [6 + t?t + % + cr6>) - (1 + 7)57?] 1 
1 + as + cr(7s — r) (1 + at) [1 + crs + a{"fs — r)] | 



2 7 + err 6* + ?/r + s(r] + aff) — (1 + 7)s?/ £ — s 

1 + as + a(js — r) 1 + crs + er(7s — r) 1 + at 



One can use equations (|2.2[) - (|2.6[) with initial values (|2.8p to derive the recur- 
rences for the polynomials Q n from several previously studied cases from [3], [4] 
and [6]. 

Here we are interested in the free harness case 7 + err = 0. After a calculation 
we get 

Proposition 2.1. Suppose a,r > 0, ar < 1, 1 + a/3 > 0, (see (|1.8p ) and 7 = 
—err. Then recurrences (|2.2|) - (|2.6|) wii/i initial condition (|2.8|) /icwe a solution which 
defines the following three-step recurrence for polynomials (Q n ) in variable y; here 
s > and x £ R are parameters. 



(2.9) 

) 1 Ci;: .r. t. si = C1 4- a^Qofn: x. t. si 4- I 

1 + crs 1 + crs 



4 crx /3 — ^(^7 crx) 
yQi(y;x,t,s) = (1 + at)Q 2 (y, x, t, s) + ( — iH — ) Qi (y; x, t, s) 



(t - s)(l + nx + ax 2 ) „ 

— Q (y; x, i, s 

l + as 
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(2.10) yQ 2 (y; x, t, s) = (1 + <rt)Q 3 (y, x, t, s) + ( a + °P) t + ft + — Q 2 ( y; ^ 4, s ) 

1 — (JT 

(t + r)(l + a/?) 

+ 7T— — TTj -Qi(y;x,t,s) 

(1 + crsj(l — (JT) 

(2.11) yQ n (y;x,t,s) = (1 + at)Q n+1 (y; x, t, s) + ( a + a ft)* + ft + aT Q n ( y; x , t, s) 

i — or 

(t + r)(l + 0!/31 . 

H t: ^ — n > 3, 

(1 — crr) z 

wf/j Qo = 1 and Qi(y; x,t,s) = y — x. 

Since p ra (y; t) = Qn(y; 0, t, 0), it is not surprising that the above recurrence coin- 
cides in this case with [31 Proposition 4.3] which we cite here for ease of reference 
in the proofs below. 

Proposition 2.2 ([3, Proposition 4.3]). Suppose (X t ) t is a quadratic harness with 
parameters such that a, r > 0, <tt < 1, 1 + a(3 > 0, and 7 = —err. If for t > £/ie 
random variable Xt has all moments and infinite support, then it has orthogonal 
martingale polynomials (p n )n given by the three step recurrences 

VPi(y; *) = (! + o-t)p 2 (y; t) + (at + /3)px(y; t) + tp (y; t), 

. . . , , , (a + o-(i)t + f3 + ar . , (i + r) (1 + afl) , , 

1 — <JT 1 — OT 

ypn{y,t) = (1 + at)p n +i(y; t) -\ Pn{y;t) 



1 — 0~T 

(t + r)(l +af3) 

(1 - (JT) 2 



Pn-i(v,i), n>3, 



with po = 1 and pi (y; t) = y. 



Thus the recurrences for the polynomials (Q n ) n and (j> n ) n turn out to be some 
finite perturbations of the constant coefficient recurrence. Therefore (Q n ) n and 
(Pn)n turn out to be Bernstein-Szego polynomials (see \12\ §2.6]), which are or- 
thogonal with respect to the probability measure with the absolutely continuous 
part of the form v 'ax 2 + bx + c/p(x), where p is a polynomial. 

3. One dimensional distributions 

Let 7T* denote the orthogonality measure of the polynomials (p n (y]t)) n with 
t > 0, given by the three step recurrences in Proposition 12.21 The existence of 
7r t is assured by Favard's Theorem. In order to examine 7r t we will compute its 
Cauchy-Stieltjes transform Gt- Recall that the Cauchy-Stieltjes transform of a 
probability measure p is an analytic mapping of the upper complex half-plane C+ 
into the lower half-plane C_ defined as 

G(z) = [ — p(dx). 

J R z — X 

One of the properties of the Cauchy-Stieltjes transform that we shall use here is 
(3.1) lim iyG{iy) = 1 

2/->oo 

(see e.g. [5]). 
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Lemma 3.1. If n 2 > 4cr > 0, 9 2 > At > and a + a/3 > then 



HI 5F • 5 Jj=°- 

Proof. It is known (see [2], (2.10) and Theorem 2.4) that with orthogonal polyno- 
mials one can associate a continued fraction (built from the coefficients of the three 
term recurrence), which, if it converges, is equal to the Cauchy-Stieltjes transform 
of the orthogonality measure of the polynomials. In the case of the recurrence from 
Proposition ^. 21 one gets 



(3.2) G t {z) = - 



where g t has the continued fraction expansion 

l 



9t{z) 



(l+«/3)(t+r) : 

1 (a+af3)t+(l+aT (l+crt)(l-<r-r) 2 

l+at Z (l+o*)(l-<rr) (l+aff)(t+r) 

1 _ (a+al3)t+l3+aT _ (l+fft)(l-gT)^ 



l+crt^ (l+ert)(l-o-r) 



1 + al 



so g t itself is the Cauchy-Stieltjes transform of a measure and it satisfies the qua- 
dratic equation 

9t ^ ~ (a + a/3)t + (3 + ar (1 + a/3)(t + r)(l + at) 

z 1 n ^2 9t{z) 

1 — (7T (1 — (TT) Z 

(To justify convergence of the continued fraction expansion of g t , one can use e.g. 
Theorem 2.1 g].) Hence 

(3.3) g t (z) = — — — 1 aT — - [ (1 - (jt)z - (a + a(3)t - (3 - ar 
2(1 + ap)(t + t)(1 + at) \ 



± 



\j ((1 - <jt)z - (a + a(3)t - (3 - arj 2 - 4(1 + a/3)(t + r)(l + at) 
By the square root in (|3 . 3[) we understand 



b + 2^/7: / 6-2^ 



az — b — 2\fc ■ y az — b + 2^ = a \ / : - : — • \ / : - - 

where a = l-<rr>0, b= (a + a/3)t + (3 + err el,c = (1 + a(3)(t + t)(1 + crt) > 0, 
the mappings z M> \J z ~ (b ± 2y/c)/a are analytic on C\{(& ± 2-y/c) /a — t : t > 0}, 
and take positive values for M 9 z > (6 + 2 v / c)/a (so they are the principal branches 
of the square root, composed with linear transformations z i-> z — (b ± 2y / c)/a). 
Hence, by (|3.1j) . one has to choose the "— " sign in the "±" in (|3.3j) . 
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Now, after inserting (|3.3[) into (|3.2p . a calculation that uses (|1.8|) gives 
rz + 6»t t[{l + (TT + 2<jt)z + tri-6] 



(3.4) G t (z) 



" rz 2 + 9tz + i 2 2((tz 2 + r]Z + l)(rz 2 + 6»<z + t 2 ) 

- <tt)z - (a + a/3)t - f3 - ar} 2 - 4(1 + at)(t + t)(1 + a/3) 



2(crz 2 + r/z + 1)(tz 2 + 6tz + t 2 ) 

Sticltjes-Perron inversion formula (see e.g. [5J Theorem 2.5 and Section 2.3]) states 
that a finite Borel measure v with the Cauchy-Sticltjcs transform G is absolutely 
continuous with respect to Lebesgue measure on the set 

A = {x : limG(a; + ie) = &(x), a finite number with Im<I>(a;) =^ 0}. 

The atoms can only be located at simple poles of G (see [2]). A very useful result 
(see [11], Chapter XIII. 6) states that if 

B = {x : lim G(x + ie) = ool, 

ej.0 

then v(M. \ (A U B)) = and v restricted to B is singular relative to Lebesgue 
measure. Therefore we see that the absolutely continuous part of 7r t is concentrated 
on the interval [a_(t), a+(i)] with 

. . (a + of3)t + (3 + ra± 2y/(l + at)(t + t)(1 + a/3) 

a ±(t) = : , 

1 — (TT 

the atoms can be located at (at most) four points being zeros of the polynomial 
(crz 2 + rjz + 1)(tz 2 + 9tz + t 2 ): 



-n ± Jr) 2 - 4cr , , 6 ± VO 2 - 4r 

c±(t) = -t- 



2a w 2r 

and 7r t does not have continuous singular part. 

Hence we will have established the lemma if we prove the following claims. 

Claim 3.2. The continuous part ofirt does not assign any probability to the interval 
[&_,&+]. 

Proof. It suffices to check that a_ (t) > b + for t > 0. It is easy to see that 



„ , . _ (l-ar)yi+^ 
[> 2{l +atf/ 2 {t + Tf/ 2 

Note that a — a/3 = r\ so 

^ (a + a(3)(l - ar) 1 + ar 
2a v / (a - a/3) 2 - Aa 2ct 
is well defined. The following bound shows that t* is in the domain of a_(t): 

(l- aT) 2 (l + a(3) 
rf - Aa 

Since a'_(£») = 0, and a_ is convex, this is a global minimum of a_. It follows that 
a-(t) > a-{t*) for t > 0. Since a_(t*) = 6+, the proof is complete. □ 



(1 + aU)(U +r)= V ^ U> 2 > ^ > °- 



Claim 3.3. The discrete part of nt does not assign any probability to the interval 
(b-,b+). 
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Proof. If 776* < then the points b± and c±(t) are separated by the interval 
(a_(i),a+(t)), so c±(t) (&_, 6+) for all t > 0. 

Suppose then that 77 and are of the same sign. Since by (|1.8|) 

77(1 + or) + 2ad 
<a + af3 = — , 

1 — 0~T 

it follows that 77 > and 8 > 0. The weights p±(i) of the points c±(t) are given by 
the residues of the Cauchy-Stieltjes transform Gt at the points c±{t) (see [5]). A 
lengthy calculation reveals that 



P-(t) 



P+(t) = 



2r -t [2r]T + (1 + (tt)(6 - y/9 2 - 4r)] + 2tV0 2 - 4t 



g+ye 2 -4T \ / . _ e+\/e 2 -4r 



<tV6> 2 - 4r (8 - V0 2 - if) U - . ,. 

\ r}—\/r} 2 —4a J \ rf-\-y/ r} 2 —Act 

-2r(-t [2r]T + (I + <tt)(6 + V¥~ If)] - 2tV(9 2 - 4t 



ctVW^P (g + V^4f) 2 ft- fl-vgEH ^ ft - e-v^Egi ^ 



where (a)+ = (a + |o|)/2. Clearly = 0, and p_(t) > only on the finite 

interval 



2ry/d 2 - 4r 

< t < 



[6 - sjQ 2 - 4f)(l + err) + 2t7t 
in particular, 



2rV6> 2 - 4t (9 + y/0 2 - 4r 

< 



(0-\/6 ,2 -4r)(l + CTT)+2?7T 77+ V 7 ? 2 - 4cr ' 

Since c_ evaluated at the right hand side of the above inequality is equal to 6+ , we 
get that the support of the discrete measure P-(i)£ c _(t) stays above the level b + 
for all t > 0. □ 

The proof of Lemma |3. H is complete. □ 

In the next lemma we briefly describe the extreme case of r = 0. 

Lemma 3.4. The assertion of Lemma \3.1\ holds when rj 2 > 4a > 0, t — 0, 9 2 > 0, 
and a + oft > 0. 

Proof. The proof of Claim I3~2l carries over to the case r = without any changes. 
Next we consider the atomic part of the measure. Instead of two lines c±, we have 
one 

t 

to take care of when examining the discrete part of 7r t . Since ?7 2 > and 2 > 0, 
as in the proof of Claim [3TB1 it suffices to consider the case of 77 > and 9 > 0. The 
residue of Gt at c(t) is 



(-« [t(i + r,e)-e 2 ]) + 

P ^ = 2o-6 2 t(t/9 + b-)(t/6 + b + )' 
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It follows that p(t) > for 

f) 2 

< t < 



l + ?70' 
in particular 



9 2 

< -Bb A 



so the support of the discrete measure p{t)S c M stays above the level b + for all 
t> 0. □ 



4. Generating functions and connection coefficients 

Our next task is to establish an algebraic relation between the polynomials (p n ) n 
and (Q n )n- In our setting, the relation takes the same form as in [SJ Proposition 
2.2]; a more complicated example occurs in [?J Theorem 2.1]. 

Proposition 4.1. There exist polynomials (bk(x,s))k and (ck(x, s))k in variable x 
such that bo(x, s) = 1 and 

n 

(4.1) Q„(y;x,t,s) = c n {x, s) + ^b n - k (x, s)p k {y;t) 

k=l 

and (bk)k an d (ck)k do not depend on t and y. 

Proof. Let Q denote the generating function of the polynomials {Q n )n, that is, let 

oo 

Q(z, y, x, t,s) = ^T z n Q n {y; x, t, s). 

71=0 

From (|4.ip with y = and t = we see that we must have c n {x, s) = Q n (0; x, 0, s). 
Therefore, to prove the proposition, we need only to verify that the right hand side 
of 

£ Q(z,y,x,t,s) - Q(z,0,x,0,s) 

Q(z,y,0,t,0)-1 
does not depend on variables y and t. Then the series expansion 

oo 

&(z, x, s) = ^z n b n (x,s) 

n=0 

defines the appropriate sequence {b n (x, s)) n >o- 

To prove (|4.2[) we need an explicit formula for Q. Using the three step recurrence 
for (Q n )n, after a routine calculation one can verify that 

(4.3) 

yQ = yQ + yzQ, + z 2 g± r) . ( | + aP) Q, + (q-Q - zQ, - z 2 Q: 

(1 + crs){l — <jt) z V 

(a + a0)t + (3 + ar /- _ \ (i + r)(l + /~ 

+ v ^ p Q-Q0-2Q1 + z , n — ^ \Q~Qo-zQ_ 

1 — (TT V / (1 — UtY V 
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(to save space, we dropped the arguments (z, y, x, t, s) in Q and (y; x, t, s) in Q n )- 
From this equation, after an elementary, but lengthy algebra, we first obtain a 
formula for Q as a rational function of z, and then verify that (I4.2[) holds true with 

b(z, x, s) = 



/ z 2 (l + aB)s + z(l - err) [s(a + a/3) - x(l - err)] + as (I - err) 2 \ 
y z 2 r(l + off) + z{B + ar)(l - or) + (1 - err) 2 + J 

Since 6(0, x, s) = 1, we get bo(x, s ) — lj as claimed. □ 

5. Quadratic harness property 

In [5] we developed an operator approach, related to Lie algebra techniques, to 
the verification of the quadratic harness property. It uses a representation of the 
process under investigation through an operator X t = X + iy, where X and y are 
some operators built from some compositions of the ^-differentiation operator D q 
and the multiplication operator Z. 

Here we show how to exploit this technique to prove the quadratic harness prop- 
erty of the Markov process with martingale polynomials given in Proposition 12.21 
Let 

(5.1) Q t * s ,„(x, y) = A Mi „x 2 + 5 MiU yx + C M:M y 2 + D t ^ u x + £ M , n y + F t , s , u 

be the quadratic form in the non-commuting variables X, y (a dual of (|1.3j) ). Define 
the generating function of the polynomials (p n ) n as 

00 

Pt( z ,v) = X] ^PniWit)- 

n=0 

In the free harness case, we are going to use the O-differentiation operator D, skip- 
ping the subscript, so 

D(g)(z) = 9{Z) = g(0) and Z(g)(z) = zg(z) 
(we treat them as the linear operators on formal series g(z) in the variable z). 
Proposition 5.1. Let 7 = —err and 

T(l + af3)~^ r(a + ad) ^ 2 ar 2 (l + aB) _3_ 2 

x = D + /3ZD+ — — Z 2 D + — — Z 2 D 2 + — ^ r^Z 3 D 2 , 

1 — err 1 — err (1 — err) 2 

_, ^„ „_9 a(8 + ar) _,9„o a/3 + crr^,9„ err(l + aB) _,q„ 9 
y = Z + aZD + erZD 2 + -Z 2 D 2 + — — Z 2 D+ -. y ZD. 

1 — err 1 — err (1 — err) 2 

T7ie operator X t = X + ty satisfies 

(5.2) X 2 = Q* S!tl (X s ,X u ) VKt<u, 
£/ie quadratic form given by (|5.ip and (11.51) . Moreover, 

(5-3) ypt(z,y) = (X t Pt) (z,y). 

Proof. A long but straightforward calculation shows that X and y satisfy the dual 
version of the ^-commutation equation 

(5.4) [x, y] 7 = erx 2 + ry 2 + ryx + 9y + I. 
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By [3 Proposition 4.9], ([Ojl holds. The algebraic identity (|Q|) follows fr om the 
three step recurrences for the polynomials {p n )n given in Proposition ^. 2\ by another 
routine calculation. □ 

Proposition 5.2. If (X t )t is a Markov process such that the random variables X t 
have moments of all orders and (p n ) n are orthogonal martingale polynomials of the 
process (X t )t, then {X t )t is a quadratic harness with 7 = — err. 

Proof. Condition (jl.ljl holds true. Indeed, 

EX t =E\p 1 (X t ;t)p {X t ;t)} =0. 

For s < t, by the martingale property (|2.1I) and the first recurrence in Proposition 
Owe get 

E[X s X t ] = E [X S E [ Pl (X t ;t)\ T s ]] = E [X sPl (X s ; s)} = 

E [(1 + <rs)p 2 (X s ; s) + (as + P) Pl (X s ; s) + sp (X s ; s)} = s. 

An efficient way to verify (|1 .4[) and (|1.2|) is to use (|5 . 3[) to represent the process 
through the operator X t from Proposition 15 . 1 1 as 

X t p t {z,X t )=X t (p t {z,X t )). 

This, together with the martingale polynomial property, which for the generating 
function p t implies 

E[p t {LX t )\X s ]=p s {tX s ), 

gives for s < t < u 

E(p s ((,X s )X t p u (tX u ))=E(p s (C,X s )X t p t (t,X t )) = 

X t E (p s (C, X s )p s (£, X.)) = X t G s (C, 0, 

where 

00 

G S (C, = E (p s (C, X s )p s (t X.)) = £(CO"E (Pn(X s ; s)) 2 , 

n=0 

and X t acts on G s (£, £) as on a series in variable £. Thus we arrive at the equivalence 
of 

(5.5) E{p s (CX s )X t p u (^X u )) =a t ^ u E(p s ((,X s )X s p u (!;,X u )) 

+ b t , s , u E (ps((, X s )X u p u (£, X u )) 

and 

X t G s (C, £) = «!...„XJ.\ (C, ?) + b t<s , u X u G a (C, 0- 

The latter (and so (|5.5l0 follows from the operator identity X t = a tt s,uX s + b tyS ,u^-u, 
which is a trivial consequence of the representation X t — X+ ty, t > 0. 
Now (|5.5j) means that 

E (p n (X s ; s)X t p m (X u ; u)) 

(p n (X s ;s)X s p m (X u ; u)) + & 4 , s , JE (p n (X s ; s)X u p m (X u ; u)) 

for all m, n > 0. Since the random variables X t are bounded, polynomials are dense 
in L 2 (X S ,X U ) (see [8j Theorem 3.1.18]). Thus by the fact that (X t ) t is Markov, 
(fl~4|) follows from (pT5]) . 
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The proof of (|1.2j) is similar, since (|1.2|) is equivalent to 

(5.6) E(p a te,X s )Xfp u (Z,X u )) = A t , s , u E(p s ((,X s )X 2 p s (£,X s )) 

+ B t , s>u E (p s ((,X s )X s X u p u (£,X u )) + C t , s , u E(£ s ((,X 8 )Xlp u {£,X u fi 
+ Dt, s , u E (p s (C, X s )X s p s (£, X s ))+ E t , s>u E (MC, X s )X u p u (t, x u)) 

+ F M ,„E (p s (C,X s )p s (^X s )). 

Observe that if s < u then 

(5.7) E($ s (t,X s )X s X u p u (Z,X u )) =E@ s ({,X s )X s X u p u (Z,X u )) = 
X u E(p s ((,X s )X s p s (t,X s ))=X u X s E(p s ((,X s )p s (£,X s ))=X u X s G s ((,(;)- 

Similarly, X 2 v G s (C,0 = E (p s ((, X s )X 2 p u (£, X u )) for v e [s,u]. This and f£Q 
show that (|5.6[) follows from the operator identity (|5.2I) applied to G s (£,£) treated 
as a formal power series in variable £, proving (|1.2|) . 

□ 

6. Construction and uniqueness 
Now we are in a position to prove the main result of this paper. 

Proof of Theorem We first note that since the time inversion (tXw t ) t of the 
quadratic harness {X t )t is still a quadratic harness with parameters rj, a replaced by 
9, t (see Remark 2.1 in [3]), it does not matter whether we construct (X t )t>o, or its 
time inversion (tX 1 / t ) t> o- Secondly, we note that if rj 2 > 4a > and 8 2 > 4r > 
then it is impossible to have simultaneously a + a(3 — and (3 + olt = (recall 
(|1.8jl ). So passing to time inversion if necessary, we may assume that a + a (3 ^ 0, 
and passing to (— X t ) if necessary, we may assume, a + a/3 > 0. Similarly, observe 
that if a > 0, r = and rj 2 > 4cr then a + a/3 ^ 0. Indeed, if r = and a + a/3 = 
then r\ = —2a6; rj 2 > 4cr implies a9 2 > 1 while 1 + af3 > implies ad 2 < 1 - a 
contradiction. Hence we may assume a + aj3 > as before. Therefore, without loss 
of generality, we will consider the following list of constraints for the parameters 
for which we want to construct the quadratic harness: 

• Case 1: u, r > and rj 2 < 4a, 

• Case 2: a, r > and rf > 4cr, 9 2 > 4r, and a + af3 > 0, 

• Case 3: a > 0, r = 0, rj 2 < 4a, 

• Case 4: a > 0, r = 0, rj 2 > 4a, and a + a/3 > 0, 

• Case 5: a > 0, r = 9 = 0, 

• Case 6: a = t = 0. 

We omit Cases 5 and 6, as the full construction of the quadratic harness with a = 
t = appeared in 6 and the case r = 9 = is the time-inversion of [5j Theorem 
4.3]. In the remaining cases, we will use polynomials (p„)„ = (p n (y,t)) n from 
Proposition ^. 21 to determine measures 7r t which will be the univariate laws of (X t ). 
The orthogonality measures P Si t(x, dy) of the polynomials (Q n ) n — {Qn(y, x, t, s)) n 
from Proposition ^. II will be the transition probabilities of (X t ) , i.e. the conditional 
laws C(X t \X s — x). We will verify that these probabilities satisfy the Chapman- 
Kolmogorov equation, so that (X t )t is indeed a well defined Markov process. 

It is clear that the coefficients at Q\ in (|2.10[) and Q n -i at (|2.11[) are nonnegative. 
So by Favard's theorem, in order to define probability measure P Sj t(x, dy), we only 
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need to check that the coefficient at Qo in ()2.9[) is nonnegative for x from the 
support of the measure 7r s . 

The coefficient at Qo is obviously nonnegative in Cases 1 and 3, as ry 2 < 4cr. By 
Lemma [3~T1 the coefficient at Qo in (|2.9[) is nonnegative in Case 2. By Lemma [3T4l 
the coefficient at Qq in (|2.9I) is nonnegative in Case 4. 

Thus, in each case, the polynomials (Q n (jJ\ x,t, s)) n determine the probability 
measures P s .t(x,dy) for all x € supp7r s . Observe that both families of measures 
(7r t ) t and (P s ,t(x,dy)) s j,x are compactly supported and uniquely determined, as 
the coefficients of the three step recurrences (|2.9j) - (|2.1ip arc bounded in n. 

We now verify that the probability measures (P s j(x, dy)) are the transition 
probabilities of a Markov process. To do so, notice that (|4.1|) for n > 1 implies 



(6.1) Q n (y;x,t, s) = ^6„_ fc (x, s) \p k (y;t) -p k (x;s)] Vx,yeR. 

k=l 

(Observe that Q n (x; x, s, s) = as a consequence of (|2.9[) - (|2.11[) .) Since b = 1 and 
p a = 1, a recursive use of (|6.1[) yields 



(6.2) p n (y;t)P Stt (x,dy) =p n (x;s) V a; e supp 7r s . 

JR 



Let 



U = < 



. \ (— J] -1 , oo), when cr = and 77 < 0, 

. \ (— 00, i]^ 1 ), when a = and 7/ > 0, 



L, when cr = = 7/. 

We proceed to show that for < s < t < u and for a set of x of 7r s -measure one 

(6.3) P s , u (x,-)= [ P t . u (y,-)P s . t (x,dy). 

Ju 

First, consider the special case s = x = of (|6.3p . which we state equivalently as 

(6.4) 7T„(-)= f P t>u (y,-)n(dy). 

Ju 

Define v{A) = f v Pt,u{Vt ^) 7r t(dy). To prove that v{dz) = -K u {dz), we only need to 
show that the polynomials Q n (z; 0, u, 0) = p n (z; u) are orthogonal with respect to 
is(dz). Since the argument is analogous to the one developed in the general case 
below, we omit it. 

From the fact that (16.31) holds for s = x = 0, we deduce that 



(6.5) P St t(x, U) — 1 Vi€ supp7r s . 

Indeed, observe first that since the coefficients in the three step recurrences (I2.9[l - 
(|2.11[) depend continuously on x, the same is true for the Cauchy-Stieltjes trans- 
forms of measures P Sj t(x, dy), which take form (13.4[) with parameters that depend 
on s,x, see (|2.8p . So U 3 x <— > P Sj t(x,U) is a continuous function. Then, Lemma 
Oand (HU) imply that 



l = 7r t (Z7) = J P s ,t{x,U)n s (dx). 

Therefore P St t(x,U) = 1 on a set of x of 7r s -probability one. By continuity of 
P s .t(x, U) in x, the conclusion follows for all x € supp7r s . 



14 



W. BRYC, W. MATYSIAK, AND J. WESOLOWSKI 



We now prove that (16. 3|) holds in general. Fix s > and x £ supp7r s , and let 
K') = lu P*.«(j/> ')^s,t(x, dy). We will show that v{dz) — P s<u (x,dz) by checking 
that the polynomials (Q n (z; x,u, s)) n are orthogonal with respect to v(dz). Since 
v is a probability measure and (Q n ) n satisfy a three step recurrence with bounded 
coefficients, to verify that v(dz) coincides with P s j(x,dz), it suffices to prove that 
v(dz) integrates Q n (z; x,u, s) to zero when n > 1. Using consecutively (|6.1I) . (|6.2I) . 
(|6.5p . again (|6.1|) . and the fact that LQ n (y;x,t,s)P s ,t(x,dy) = Q for n > 1, we 
get 



Thus (|6.3p holds and P Sj t(ir,dy) are transition probabilities of a Markov process 
(X t )t with state space U. Since p n {y',t) — Q n (y;0,t,0) it follows from the con- 
struction that for fixed t > polynomials (p n (y',t)) n are orthogonal with respect 
to 7Tt(dy) = Po,t(0,dy); their martingale polynomial property follows from (|6.2|) . 
Proposition 15.21 implies that (Xt)t is a quadratic harness with parameters rj, 9, a, 
t and 7, with 7 = —err. 

Uniqueness of the process (X t ) t follows from the fact that orthogonal martingale 
polynomials (p n ) n determine uniquely the joint moments of the process. Recall 
that the measures Tr t are compactly supported, so the joint moments determine the 
finite dimensional distributions of the process uniquely. □ 
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